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Abstract. We review several techniques and ideas initiated by a remarkable 
work by Spencer 1261 . used and further developed in numerous subsequent 
researches. We also describe a relatively short and elementary derivation of 
the spectral and strong dynamical Anderson localization from the fixed-energy 
analysis of the Green functions, obtained either by the Multi-Scale Analysis 
(MSA) or by the Fractional-Moment Method (FMM). This derivation goes in 
the same direction as the Simon- Wolf criterion 1281 . but provides quantitative 
estimates, applies also to multi-particle models and, combined with a simplified 
variant of the Germinet-Klein argument 1201 , results in an elementary proof of 
dynamical localization. 



1. Introduction. 

The mathematical theory of Anderson localization, describing the motion a quan- 
tum particle (or a collection of non-interacting particles) in a disordered environ- 
ment has reached by now its age of maturity. The number of different mathematical 
models (different even from the point of view of applications to physical systems) 
and technical tools, allowing to analyze these models, is quite impressive. On the 
other hand, this also makes mastering these techniques difficult for the beginners. 
It is not always easy to see the main ideas behind dozens of pages filled with defi- 
nitions, preliminary facts and sophisticated arguments pushed to the extreme due 
to the complexity of the problem. Yet, the maturity of a mathematical theory 
can also be judged by the presence of comprehensive techniques and simple general 
principles, guiding one through the jungle of more complex models and methods, 
so a number of leading researchers in this area of mathematical physics have been 
conducting a quest for simpler, more elementary proofs of localization, intuitive (yet 
rigorous) techniques and principles. 

The quest began already in late 1980's, when Simon and Wolf [28] proved that 
suitable fixed-energy bounds on the Green functions imply a.s. pure point spec- 
trum, and then Spencer proposed in a remarkably short paper |26j an elementary 
reformulation of the fixed-energy MSA developed in his pioneering joint work with 
Frohlich [TB] . Besides the fact that the paper [55] came as a perfect complement to 
the Simon- Wolf argument, it draws the reader's attention to the parallels between 
the theory of random operators and a more traditional spectral analysis of almost 
periodic operators, following the classical works on periodic operators. 

The first paper on the FMM, published by Aizenman and Molchanov [T] in 1993, 
a few years after the cycle [THESITTJISHIHE] i a l so h& s been an important event in 
the quest for elementary proofs (which the title of jT] makes explicit). It took some 
time and efforts to complete the first stage of the FMM (a fixed-energy analysis 
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of fractional moments of the resolvents) with additional arguments leading - in a 
simple way - to the strong dynamical localization. One can only regret that the 
MSA and the FMM have been evolving in parallel with a very weak interaction, 
over the last two decades. 

In the present paper, I am going to give a short and certainly incomplete review 
of techniques and ideas brought to life or initiated indirectly by Tom Spencer in his 
short article [26] which has the good fortune and privilege to belong to those works 
which " tell much more than they say" . 

The present text is intentionally left relatively short, to keep up with the spirit of 
[21)] (which was 10 page short). For this reason, interesting applications of the MSA 
to DSO with quasi-periodic (and more generally, deterministic) potentials had to 
be omitted; these include deep analytic works by Bourgain, Goldstein and Schlag 
(cf., e.g., [1113), where the MSA techniques were applied to lattice Schrodinger 
operators with "analytic" potentials, and a recent article [7J where a different MSA- 
based approach has been used to treat parametric families of almost periodic and 
some other deterministic operators by traditional - and very simple - methods of 
the theory of random operators. On a personal note, I have to acknowledge with 
pleasure numerous fruitful discussions on this subject with Tom Spencer. 

The pioneering work [17] on Anderson localization in multidimensional disor- 
dered media (following the work by Frohlich and Spencer [16] ) proved only the 
exponential spectral localization, i.e., pure point spectrum and exponential decay 
of generalized eigenfunctions. The relations between the spectral and dynamical 
manifestations of the Anderson localization phenomenon have been studied later 
in the work [24] which has influenced further development by Germinet-De Bievre 
[T§] and by Damanik-Stollmann [13 . The overall result of these researches (per- 
fectly summarized by the title of the paper [13] ) was a clear understanding that the 
variable-energy MSA (VEMSA, in short) provides a sufficient input for the proof of 
the strong dynamical localization, in discrete and continuous random media. 

Germinet and Klein [20] made a further step and gave a much shorter derivation 
of the dynamical localization from the key MSA estimates, avoiding a tedious analy- 
sis of the random geometry of the so-called centers of localization (the latter notion 
essentially goes back to [24 ). The elegance of their approach resides, in particular, 
in the fact that the eigenfunction correlator bounds are inferred from those provided 
by the MSA directly in the entire configuration space (a Euclidean space, in their 
case). However, this elegance comes with a price: one has to rely upon a deep anal- 
ysis of weighted Hilbert-Schmidt norms of spectral projections of Schrodinger (or 
some other) operators in L 2 (R d ). Such an analysis has been carried out by Simon 
[2"5] for potentials bounded from below by — C(\x\ + 1), and later extended by Po- 
erschke and Stolz [27] to potentials bounded from below by — C(|a;| 2 + 1); the latter 
is a usual condition for the essential self-adjointness of the respective Schrodinger 
operator. In arbitrary finite volumes, the analog of Germinet-Klein argument is 
reduced to a three- word instruction: " Apply Bessel 's inequality 1 ' 1 . 

The situation is particularly simple for operators on a countable graph, where 
functional analysis is in fact replaced by linear algebra. Indeed, the required eigen- 
function expansions are completely elementary here, for Hermitian operators in 
finite-dimensional Hilbert spaces. In continuous configuration spaces (Euclidean 
spaces and quantum graphs), a similar effect is achieved whenever the random 
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operators in question have (as they usually do) compact resolvent in any finite vol- 
ume, under reasonable requirements on its geometry. Physically speaking, only such 
finite- volume bounds (uniform in the size of the volume) are relevant for applications 
to the quantum transport in disordered media. The term "finite volume" should 
not be misleading: a sample of a random media of the size of the Milky Way is still 
finite ... and the task of designing computer processors (let alone nano-devices) of 
that size does not occupy yet the minds in the physics community. 

Mathematically speaking, once uniform bounds on eigenfunction correlators are 
obtained in finite volumes of arbitrarily large size, one is just three words away from 
the strong dynamical localization in the entire space: "Apply Fatou lemma" (cf., 
e.g., 0, ®). 

With these observations in mind, I propose here a streamlined derivation of 
the VEMSA-type probabilistic bounds from their simpler FEMSA counterparts. 
Several elements of such a derivation appeared earlier, e.g., in [53], [2] (but the 
approach from Section [7] appears to be original). It plays a role similar to that of the 
Simon-Wolf argument, providing a 'soft way' from the fixed-energy localization to 
stronger manifestations of the Anderson localization phenomenon. The key notion 
here is what is called in Sections HH7] the "singular width" of the spectrum: the 
total Lebesgue measure of a (reasonably large) finite number of intervals chosen 
in such a way that outside this "singular zone" Green functions are exponentially 
small. (Curiously, the abbreviation "SW" suits to " Simon- Wolf ' , "singular width" 
and " soft way" ; the central symetry transforms it into " MS" , as in " Martinelli- 
Scoppola", cf. [33] ■) 

The structure of this paper is as follows. 

• Main notions and notations are introduced in Section [5] 

• In Sections [3HIJ we give a streamlined version of the analytic core of the 
fixed-energy MSA, in the simplest form going back to [33], but formulated 
in a more abstract fashion. 

• In Section[5j we show that the simple, fixed-energy analysis from |26j can be 
easily improved so as to provide the key probabilistic bounds on the Green 
functions stronger than any power law. 

• In Section [5] following essentialljQ [H] , we derive from fixed-energy bounds 
their variable-energy counterparts. The obtained results apply also to the 
FMM-type bounds (which are always initially obtained at fixed energy). 

• In Section[71 we give another derivation of the variable-energy bounds which 
allows to prove directly the exponential spectral localization (i.e., the ex- 
ponential decay of eigenfunctions), as well as the dynamical localization, 
under a stronger assumption on the random potential (cf. Eqn (|7.2[l ). 

• In Section[8] we formulate a finite- volume variant of the Germinet-Klein ar- 
gument, allowing to derive the strong dynamical localization from variable- 
energy bounds on the resolvents (obtained by the MSA or by the FMM). 

• Section [5] describes a simple adaptation of the techniques from Section [5] 
sufficient for the proof of sub-exponential dynamical localization. 

• In Section 1 1 01 we briefly describe another simple approach (developed in 
our recent paper [8]) which also has its roots in [26] • However, the main 
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object of the scale induction is here the decay of the eigenfunctions in finite 
balls, rather than the decay of Green functions. 
The principal statements are theorems 16.11 16.21 18.11 19.51 and 19.61 
For the sake of brevity and clarity of presentation, I do not discuss several power- 
ful (but more complex) techniques from the works by Germinet-Klein, including the 
bootstrap MSA (cf. [50]) and spectral reductions from [2T] used in a very general 
framework of random operators in R d with singular probability distributions. 

To conclude the introduction, I would like to emphasize the role that the paper 
Prl] has played in a recent development of the multi-particle MSA (MPMSA). In our 
joint works with Yuri Suhov (cf., e.g., [B]), we aimed initially to prove the spectral 
localization, which requires traditionally a variable-energy analysis. However, the 
fixed-energy analysis has a substantial advantage to simplify both geometrical and 
analytical ingredients of the MPMSA. I plan to address this subject in a forthcoming 
work, using the reductions described in Sections HHE1 

2. Basic notations, facts and assumptions 

Throughout this paper, we work with discrete Schrodinger operators (DSO) act- 
ing in Hilbert spaces of square-summable complex functions on connected countable 
graphs. Indeed, the techniques and results of the MSA, initially developed for op- 
erators on periodic lattices, are naturally extended to more general graphs with 
polynomially bounded growth of balls (such graphs as Bethe lattices remain so far 
out of the MSA's reach). Another motivation for presenting the new approach on a 
graph comes from the fact that the natural language for the description of a system 
of N > 1 interacting indistinguishable quantum particles (bosons or fermions) is 
that of a symmetric power of the configuration space Z of the respective single- 
particle system; already in the case where the configuration space is Z = Z d , d > 1, 
its A-th symmetric power is no longer a periodic lattice. 

Consider a finite or countable connected graph (G,£), with the set of vertices 
G and the set of edges £ ; for brevity, we will often call Q the graph, omitting the 
reference to £ . We denote by dg(-, •) (sometimes simply by d(- , •)) the canonical 
distance on the graph Q: dg(x,y) is the length of the shortest path x ~» y over 
the edges. We will assume that the growth of balls Bl(x) := {y : dg(x,y) < L} is 
polynomially bounded: 

sup|B L (x)| < C d L d , L>1. (2.1) 
xeg 

In particular, the coordination number ng(x) := {y : dg(x, y) — 1} of any vertex 
x is bounded by Cd (even by d — 1). The canonical (negative) graph Laplacian 
(— Ag) on a finite or countable graph (G,£) is given by 

(-A g f)(x) = J2 (/(*) - f(v)) = rig(x)f(x) J2 ( 2 - 2 ) 

where we use a popular notation (x, y) for a pair of nearest neighbors x,y £ G, i.e., 
dg(x,y) = 1, and ng(x) is the coordination number of the point x. For brevity, 
we will sometimes use slightly abusive notations like (x, y) £ A, A C G instead of 
(x,y) e (A x A)n£ g . 

From this point on, unless otherwise specified, we will use the notation G only for 
finite connected graphs, while Z will stand for a countable connected graph with 
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polynomial growth of balls. In operator form, we can write 

-Ag = rig - 22 T^b, T XjV = \l x )(l y \, 

{'■c,y) 

where ng is the operator of multiplication by the function x i— > ng(x). Given a 
proper (connected) subgraph ACg, define its internal, external and the so-called 
edge boundary (relative to Q) as follows: 

d g A = {yeA: dg(x,Q\A) = 1}, dgA = dgQ \ A, 

a e A = {(x, 2 /)ea e AxO+A: d e (x,y) = l}. 

Working with a given graph g(c Z), we always mean by a ball B^(m) C Q the 
set dg(u,y) < R}, i.e., the ball relative to the metric space (Q,dg). 

The Laplacian (hence, a DSO) in a subgraph Act/ can be defined in various 
ways. The two most popular choices are: 

• The canonical (negative) Laplacian in A, (—A^f) = (— Aa/), defined as in 
(|2.2p with Q replaced by A. It this context, it is usually considered as an 
analog of the Neumann Laplacian, and reads as follows: 

(-Alf)(x)=n A (x)- J2 /(f)- ( 2 - 3 ) 

{x,y)€A 

• The Dirichlet Laplacian (—A Ag ) = 1a(— A^ e ) 1a \ £ 2 (A). Here we use a 
natural injection £ 2 (A) £ 2 (G). The Dirichlet counterpart of (|2.4[) is 

{-AP K f){x)=ng{x)- f(v)> ( 2 - 4 ) 

(x,y)eA 

with ng(x) > nA(cc), so (— A°J > (— AjJ) in the sense of quadratic forms. 
We will use the Dirichlet Laplacians and DSO H A . Given a decomposition 
G = A U A c , A c := g \ A c , we can write 

— Ag = ng — ^ T x ^y — ^ T x>y — (^x.a + ^y,x) 

{x,y)eA {x,y)eA c {x,y)GdA 

= ((-A°)e(-Aj.))-r Al0 

with L A , g = E( a) a>e9A ( r *,a + r y,*)- Respectively for the DSO Hg = -Ag + V, 
where V : g — > K is usually referred to as the potential, one has 

Hg = H' A - T A>g , H' A := (-Ag + V) 8 (-A& + V). 

We omit the superscript " N" , since the nature of the boundary conditions in g is 
not related to the choice of Dirichlet or Neumann decoupling induced by g = AU A c . 

The spectrum of a (finite-dimensional) operator Hg, i.e., the set of its eigenvalues 
(EVs) counted with multiplicities, will be denoted by E(iJg). 

In a number of formulae and statements, we will use the parameters /3,t, g G 
(0, 1), and a £ (1, 2). Unless otherwise specified, we assume that (3 = 1/2, r = 1/8, 
g = (a— 1)/2 = 1/6 and a = 3/2. Note that the exponent figuring in Definition 
12. II then equals 7/8. (These settings will be changed in Section [5J) 

Definition 2.1. Given numbers £el and m > 0, a ball B^(u) is called 

• _E-resonant (E-K, in short), if dist(E(i?Bi(«)), E) < e~ hP , and i?-nonresonant 
(.E-NR), otherwise; 
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• (E, m)-nonsingular ((E, m)-NS), if for all x,y £ Bi(w) with d(x,y) > 

C d 2 L d - |G Bji(t0 (*,!/; £)l < e^™^^, (2.5) 

where 

7(m,L) := to(1 + L~ r ), (2.6) 
and (_E, m)-nonsingular ((£',m)-NS), otherwise. 
Observe that for any ball B L (u), \dB L (u)\ < CjL d , by virtue of ([2T]i . 

2.1. Assumptions on the random potential. For clarity of presentation, we 
always assume that the random potential field l/:ZxSl— s-lRona graph Z is IID, 
with Lipshitz continuous marginal probability distribution function (PDF) Fy : 

sup(F v (t + s) - F v {t)) < C w s, C w € (0, +oo). (2.7) 

tea 

It is well-known that this assumption can often be relaxed to uniform Holder con- 
tinuity, and even to a form of log-Holder continuity. 

2.2. The Wegner estimate. The original result by Wegner [33] has been adapted 
to a large number of classes of random operators. Here we apply its simplest version, 
for DSO with a Lipshitz continuous IID random potential. The proof can be found 
in a number of books and review articles; cf., e.g., Lemma VIII. 1.8 in |10j . 

Lemma 2.1 (Wegner estimate). Under the assumption (12.7[) . for any finite graph 
Q of cardinality \Q\ and any e e [0, 1] 

supP{dist(E(F s ),£J) < e} < C w \G\e. (2.8) 

In fact, the above statement remains valid for any ensemble of random operators 
of the form V( - ;uS) + Hq with a non-random operator Hq, for only the diagonal 
part V : Q x fl — >• K is used in the proof (cf. [10]). 

3. Decoupling of resolvents on graphs 

3.1. Geometric resolvent inequality. The second resolvent identity implies the 
so-called Geometric resolvent equation for the resolvents Gg(E) = (Hg — E)^ 1 , 
G A c(E) = (H A „ - E)-\ G gjA (E) = (H' A - E)~h 

Gg(E) = G'g(E) + Gg(E)T A ,g Gg(E). (3.1) 

For x,u € A and y € A c , one has Gg(x,u;E) = G^(x,u; E) and Gg(x,y;E) = 0. 
This results in the Geometric resolvent equation for the Green functions 

Gg(x,y;E)= £ G%(x, u; E) Gg(u\ y; E) (3.2) 

{u,u')edGA 

and the Geometric resolvent inequality (GRI) 

\Gg(x,y;E)\< ]T \G%(x, u; E)\ \Gg(u', y; E)\. (3.3) 

<u,u')e9SA 
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4. SUBHARMONICITY ON GRAPHS 

4.1. Regular subharmonic functions. 

Definition 4.1. Let Q be a finite connected graph, L > £ > two integers and 
q G (0, 1). A function / : Q — > R + is called (£, <7)-subharmonic in a ball B/,(it) C 
if for any ball B^(x) C Bl(u) one has 

/(x)<g max f{y). (4.1) 

We will often use the notation A4(f, A) := max l£ A 1/0*01- 

Lemma 4.1. //a function f : Q — > R + defined on a finite connected graph Q is 
(£,q)- subharmonic in a ball 3l(%) £ Q, with L > £ > 0, then 

f(x)<ql^\M(f,G)<q^M(f,g). (4.2) 
In fact, the factor M(f, Q) in the RHS of (|4.2I) can be replaced by M(f, B L+1 (x)). 



L 

C+ 



j\ > 1- Set A, := B j(e+1) (x), 



Proof. Since L > £ > 0, we have n + 1 : 

< j < n, and note that A n+ i C Bl+i{x), since (n + 1){£ + 1) < L + 1. Further, 
if y G Aj with < J < n and z € B£ + i(y), then, by triangle inequality, z G A J+ i. 
Consider a monotone non-decreasing function h : [0,L + 1] n N n- K + defined by 
/i(r) = M(f,B r (x)). Using the (£, g)-subharmonicity of /, we obtain 

h{j{£ + 1)) < gmax max /(^) < qh((j + + 1)), 

yeAj zeB i+1 ( v ) 

in particular, 

h(n(£ + 1)) < + 1) < qM{f, Q). 

Since /(x) = /i(0), the claim follows by the backward recursion in j from n to 0, in 
n steps. □ 

Example. Q = [0, R] D Z, R = L + 1, £ = 0, and f : x ^ q R ~ x , x E Q. For all 
y G B L (0) = [0,L], one has 

/(.t) = g^ 21 = q max <j fl ~ y = q max /(y), 

|il-x|<l |i/-a:|<l 

which implies the (0, g)-subharmonicity of the function / in B^(0). In fact, here 
the inequality of the form (|4.1D turns out to be an equality, and one has 

f(0) = q L+1 = q^f(L + l), 



which shows that the estimate from Lemma 14.11 is sharp. Note also that the in- 
equality (|4.1[) cannot be extended to the exterior point y = L + 1, since 

f(L + l) = l>q = f(L). 

Clearly, a function (£, (7)-subharmonic everywhere in Q, with q < 1, must be zero: 

< max f(x) < gmax/(j/). 

x y 



Naturally, the main raison d'etre of the Definition 14.11 is the following fact 
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Lemma 4.2. Consider a finite connected graph Q and a ball B^(u) C Q 7 with 
L > £ > 0. Fix numbers E G R, m > and suppose that all balls Be(x) inside 
Bi(ti) are (E, m)-NS. Then W y £ G \ Bi(w) the function 

f:x^\G g (x,y;E)\ 

is (£, q)-subharmonic in Hl(u) with q = e~i( m ' e ) e _ 

Proof. The claim follows directly from the Definitiori l4.il □ 

Lemma |4~T1 suffices to assess the Green functions in a ball Bi(u) which does not 
contain any singular £-ball, but to analyze the situation where Bi(it) contains one 
singular ball Bi(w) (more precisely, it does not contain any pair of disjoint singular 
£-balls), one needs the following extension of Lemma [4.11 which exploits the idea 
used of the proof of Theorem 1 in approaching a single "bad" ball separately 
from the points x and y. 

Lemma 4.3. Let Q be a finite connected graph, and f : Q x Q — > M + , / : (x,y) \— > 
f(x,y), be a function which is separately (£,q)-subharmonic in x £ B r '(u') C Q and 
iny S B r » (it") C Q, with r', r" > I > and d(u', u") >r' + r" + 2. Then 

/(«',«") < aM + ^J.M(/,g x G) < q^^Mif^g x g). (4.3) 

Proof. For each y" € B r // + i(u") define the function f y n : x' M- f(x',y") in g. By 
assumption, it is (£, g)-subharmonic in B r /(u'), so Lemma \4 . 1 1 implies . 

Vy" e B r „ +1 (u") f(u',y") = /„«(«') < q^M(f,gx Q). 

Consider now another function, /„» : y" t— > /(u', y"), y" G 5. It is {£, g)-subharmonic 
in B,,"(m"), by hypothesis. The above inequality reads as 

M{f u ,,B r „ +1 [u")) < q ^M(f,g x g), 

so another application of Lemma 14.11 proves the claim: 

/(«',«") = /„/(«") < ?wM(/„sB rni (/)) < ?TM(/,g x g). □ 

5. Fixed-energy scale induction 

5.1. Scaling of Green functions in absence of tunneling. 

Definition 5.1. A ball B£ t+1 (it) is called £?-tunneling (E-T) if it contains two 
disjoint (E, m)-S balls of radius Lk, and iS-non-tunneling (_E-NT), otherwise. 

Lemma 5.1. // a ball B Lk+1 (u) is E-NR and E-NT, then it is (E,m)-NS. 

Proof. (See Fig. 1.) Fix two points x,y € Bi t+1 (i/) with d(x,y) > L^, +£l = L^ 6 
and let R = d(x,y), so R - 2L k > R(l - 2L^ e ). Since B Lfc+1 (u) is E-NT, there 
is a ball B2L fc (w) such that any L^-ball disjoint with B2L fc (w) is (£7, m)-NS. By 
triangle inequality, there are integers r',r" > such that r' + r" > R — 2Lk, the 
balls H r i(x) and B r "(y) are disjoint from each other and from B2L (e (u>), so any ball 
Bi fe (u) inside B r >(x) and inside B r »(y) is (£7, m)-NS. 

Assume first that r' > Lk and r" > Lk (otherwise, one of the points x, y is 
covered by B2L fc (iu), so one of the radii r' ,r" > R — 3Lk and the same argument 
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B r „(y) 




Figure 1 . Example for the proof of Lemma 15.11 

as below applies). By Lemma [4.21 the function / : (v,z) n- \Gb l (v, z; E)\ is 

(Lk, g)-subharmonic in v € B r i(x) and in z g B r »(y), with q < e- mR ( 1 + L k ). By 
Lemma |4.3[ one can write, with the convention — InO = +oo: 

-kif(x,y)>-ln{ ( e -Mi+LD^) L " +1 e L *\ 



> mR ((1 + §ir) ifti C 1 - 2L D - ) 



(with m > 1, i? > i^. +e = d(x, y) > L 7 / 6 = 1/2 and Lq large enough) 



> mR | 



((l + L-^^-SL-^-L- 2 / 3 ) 

>md(x )2 /)(l + iL- 1/8 ) 
> j(m,L k+1 )d(x,y) + ln\dB Lk+1 (u)\, 
as required for the (E, m)-NS property of the ball Bx, fc+1 (u). □ 

5.2. Scale induction. Introduce the following notations: 
P k = sup P { B Lk (x) is (E, m)-S } , 

Q fe = 2 sup P { B Lfc (x) is £-R } < 2C W C d L d k e~ L " 
xez 

(the latter inequality uses the Wegner estimate ([2 
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Theorem 5.2. If there is an integer Lq > 1 such that (|2.7p is fulfilled and 

min{P ,Q } < CfL^ K , k>^, 
then for all k>0, P k < C*~ 2 L^ K . 

Proof. (The main argument combines the ideas from |j26j and |20].l By virtue of 
Lemma T5. 11 if a ball Bi fc+1 (u) is (E,m)-S, then it is either E-R or E-T. There are 
< \C\L\f\ rl pairs of disjoint L^-balls in B£ fc+1 (ii), thus 

Pk+i < -jClLl d + {Pl + -Qk+i- 



By Wegner estimate (|2.8[) . Qfc+i < CwCdLf +1 e L *+i. An elementary calculation 
shows that the function 

/ : L i-^-ln (const - In (const L^e -1 ^ 2 ") = L 1/2 - const In L 



on [l,+oo) is cither non-negative or admits a unique zero. In either case, the 
assumption Q < Cj 2 L^ K implies Qk+i < ^4 L k +i for all k > 0. Thus by 
induction on k = 0, 1, . . ., 

n-+i S ^ • 2 G d^fe+i^fc S . 2 L k+\ - °rf ^fe+r 

provided that ^ + 2d> k, i.e., for all k > □ 

This calculation shows that the power-law bound on the probabilities Vk is re- 
produced at each scale L k with the same exponent n > 0. Actually, it comes with 
a small bonus which seems to have never been used, until recently. Pick a value 
K > l 2 ^, so that :=(--—)- 1 > 0, and observe that 

M1±^L _ 2d > MiM! - 2d = «(1 + 0) fe {£ - f } = k(1 + 

If for some fco > one has P ko < C^ 2 L k( ^ 1+e ' > ° , then by induction, for all k > ko, 

S 2 • ^dl^k+l^k S -jj4 L k + l 
2 T -K(l+8) k + 1 



- C d L k+X 

This gives rise to the probabilities P k decaying faster that any power law. 

Observe also that taking ail results in 2a/(2 — a) 1 2, so that the exponent 
k > in the hypothesis of Lemma 15.21 can be chosen arbitrarily close to 2d. 

We see that the complete, fixed-energy MSA procedure can be effectively reduced 
to Lemma I5.il and Theorem \5.2\ and results in upper bounds on the probability of 
singular balls decaying faster than any power-law. 

6. From fixed to variable energy: First approach 

Now we establish a fairly general relation between fixed-energy probabilistic es- 
timates on the Green functions and variable-energy bounds for two disjoint finite 
volumes. It does not matter how the probabilistic input is obtained; in particular, 
the results of this section can be combined both with the MSA, performed for each 
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fixed energy E in a given interval J C K, and with the FMM (which always starts 
as a fixed-energy analysis). 

It is convenient to assume that |7| = 1, so the interval I with the Lebesgue 
measure mes(-) is a probability space, and so is the product space (fi X J,P X mes). 
(The idea of using the " disorder-energy" space with product measure has been used 
earlier, e.g., in [52] and [25].) Given LeN and points x,y € Z, set for brevity 

M x . y {E) = \G BL (x,y;E)\, M X (E) = max M x . y (E), (6.1) 

and introduce the subsets of I parameterized by a > 0: 

£ x , y {a) = {EeI: M x , y (E) > a}, £ x (a) = {E E I : M X {E) > a}. (6.2) 
(The L-dependence will be often omitted for brevity.) 

Theorem 6.1. Let L > 0, x G Z, y G 8~Bl(x). Let {Aj}jLi be the eigenvalues of 
the operator H BL r u \(u)) an d / don interval of unit length. Let be given numbers 
a, b, c, Vl > such that 

b<min{\B L (u)\- l ac 2 , c} , (6.3) 

and for all E € L 

P{£ x (a)} = P{M x (E)>a}<V L . (6.4) 
There is an event B x (b) with ¥ {B x (b) } < b Vl such that V u> £ B x (b), the set 

£ x (2a) = £ x (2a;w) = |e : M X (E) > 2a| 

is contained in a union of intervals \JjLj_Ij, Ij := {E : \E — Xj\ < c}, Xj G /. 
Proof. Consider the following events parameterized by b > 0: 

B x (6) = {u>efi: mes(£ x (a)) >&}. (6.5) 
Apply Chebyshev's inequality and the Fubini theorem combined with (|6.4[) : 
P{B x (b)} < 6- 1 E[mes(£: :E (a))] 

= 6- 1 J^dEE[ l {M AE)>a} ] < b-'ViL). (6.6) 

Now fix any uj B x (b), so that mes(£ x (a); ui) < b. There is a subset {Aj};^ of 
the EVs of the operator H Bl m such that the Green function E h-> Gb^)^, y; -E) 
reads as a rational function (below we remove the vanishing terms, if any) 

N' 

f:E^G BL{x) (x >y ;E)=:J2j^rE> N' < N := \B L (x)\; (6.7) 

j=i J 

here Kj = Kj(x,y) ^ and £V < £V |V>i(£)-*M2/)l < Let 
7£(2c) = {AeR: min |Aj - A| > 2c}, 

K(c) = {AeR: min |Aj - A| > c}, c> 0. 

Observe that, with < b < c, A b := {E : dmt(E,K(2c)) < b} C 7£(c), hence, the 
set Ab is a union of open sub-intervals at distance > c from the spectrum, and on 
each sub-interval one has \ f'(E)\ < Nc~ 2 . Let us show by contraposition that, with 
u) B XlV (b), 

{E: \G BL {x,y;E)\>2a}nK{c) = 0. 
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Assume otherwise, pick any point A* in the non-empty set in the LHS, and let 
J := {E' : \E' - A*| < b} C A b C 11(c). Then for any E € J one has, by 



|/(£)| > |/(A*)| - | J\ sup |/'(£')| > 2a - ^Vc- 2 • b > a, 
E'ej 

so J C £ Xt y(a) and mes(£ Xi!/ (a)) > mes(J) = 2b > b, contrary to the choice of w. 
Since the set 1Z(c) is independent of y, the assertion follows from (|6.6[) . □ 
Below we provide some possible choices of the parameters a, 6, c (depending, of 
course, upon the scale L) in three most frequently used frameworks. 

(1) Weaker MSA-type bounds: T Lk = Pk < L k K(1+9)k , n > a = 3/2. 
One can set, for L E {L k , k > 0}, 

a{L k ) = L k 5 , b{L k ) = L k -> , c(L k ) = L k yo 

(2) Sub-exponential MSA-type bounds: "P Lfc < e~ L S 5 > 0. Then one can set 

a(L fc )=e-^, b(L k ) = e'i^ , c(L k ) = e^ L *. 

(3) FMM-type bounds: V L < e~ mL , m > 0. Then one can set, for all L € N 
large enough, 

a(£) = e-3 mi , &(L) = e-i mL , c(L) = e-t L . 

Theorem 6.2. Assume the condition (|2.7[) . Suppose that for some LgM, numbers 
a = a(L), b = b(L), c = c(L) and Vl > o&ey (|6.3[) . and /or some interval Jcl 
and aZZ E E I, for any ball B^(a;) C Z 

F{M X (E) >a}<V L - (6.8) 

T/ien /or any two disjoint balls B^(x), B^{y) c -H £/ie following bound holds true: 

F{3EeI: mm(M x (E),M v (E)) > a (L)} < AC w C 2 d L 2d c(L) + (6.9) 

Proof. Let the events B x (b), B y (b) be defined as in (|6.5[) and Z? = B x U then 

F{£ x (a) n£ y (a) ^ 0} <F{B} + E [P{ {£ x (a) l~l £ y (a) ^0}n6 c }] 

< 26"Vl + E [F { {£ x (a) n /0}nB c | J^) } ] . (O0) 

It remains to assess the conditional probability in the RHS. For oj g" B c , each of 
the sets £ x (a), £ y (a) is covered by intervals of width 2c(L) around the respective 
EVs \i(x) € S(i?B £ (x) 5 ^j(y) € S(i/ Bl ,(y), and for disjoint balls Bl(x), Bl(v) these 
spectra are independent. Now apply Theorem 16.11 and the Wegner estimate: 

P{{£ x (a)n£ y {a)^0}nB c \F BL{y) } < \B L (y)\ su P P{ dist(^(a), A) < c} 

< \B L (y)\ supP{dist(£(# B£ (s),A) <2c} 
xei 

<4C w \B L (y)\-\B L (x)\-c(L). 

(6.11) 

Collecting (|6.10p and (|6.11l) . the assertion follows. □ 

It is clear that the above approach, albeit very general and based on a Wegner- 
type bound, gives rise to exponential decay bounds on the Green functions only if 
the fixed-energy analysis provides exponential probabilistic bounds; as it is well- 
known, this has been achieved so far only in the framework of the FMM. 
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6.1. Spectral localization. The assertion of Theorem 17. 3 1 has a structure similar 
to that of the MSA bound from the work by von Dreifus and Klein [T^]. More 
precisely, it guarantees a decay rate of Green functions slower than exponential, 
but faster than any power-law. It is not difficult to adapt the well-known argument 
from |12j and prove that with probability one, all polynomially bounded solutions 
to the equation H(ui)ip = Ei\) are in fact square-summable. The latter property 
requires a Shnol-Simon type result on spectrally a.e. polynomial boundedness of 
generalized eigenfunctions. It will follow independently by RAGE (Ruelle-Amrein- 
Georgescu-Enss) theorems (see a detailed discussion along with a bibliography, e.g., 
in [9]) from the dynamical localization proven in Section [H 

7. From fixed to variable energy: Second approach 
7.1. The spectral reduction. 

Theorem 7.1. Let be given a ball 3l(x), L > 1, and numbers a(L), b{L), c(L), 
Vl > obeying (|6.3[) and such that, for some interval I . all E G I , 



P{M x (E)>a} <V L - (7.1) 
Then the following properties hold true: 

(A) For any b > Vl there exists an event Bb such that P{Sb} < b~ x Vh and for 
any uj G" B v the set of energies 

£ x (a) = £ x {a; u) = {M X (E) > a} D I 

is covered by K < 3N 2 intervals Ji = [Ej~ , Ef] of total length | J{\ < b. 

(B) The endpoints Ef 1 are determined by the function E t— > (1^. | [H-q l ^—E)^ 1 \ l y ) 
in such a way that, for the one-parameter family A{t) :— H^ L t u \ +t 1, the end- 
points E i (t) for the operators A(t) (replacing operators Hb l (u)) have the form 

E±(t) = Ef(0)+t, tew. 

Proof. (A) Fix a point y G d^B^^x) and consider the rational function 
fy - E ^^X^E- = ^ K-E ■ 

l—l 1=1 

Its derivative has the form 

and has < 2N — 2 zeros and < N poles, so f y has < 3iV — 1 intervals of mono- 
tonicity Jj^, and the total number of monotonicity intervals of all functions {f y , y G 
d~B L (x)} is bounded by K < |<9~B L (x)|(3A - 1) < \B L {x)\{ZN - 1) < 37V 2 , so 

U yed - BL{X) {E : f y (E) >a} = U^Ji, J t = [Er,E+] C I, 

where, obviously, J2i \J*\ — mes{£' : M X (E) > a}. 

(B) Consider a one-parameter operator family A(t) = H Bl / u \(ui) +tl. All these 
operators share common eigenvectors; the latter determine the coefficients K{, so 
one can choose eigenfunctions ipiit) constant in t and obtain Ki(t) = ftj(O). The 
eigenvalues of operators A(t) have the form \(t) — Aj(0) +t. We conclude that the 
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Green functions, with fixed x and y, have the form f x ^ y (E;t) = f x , y (E — £;0), so 
that the intervals Ji(t) have indeed the form J,(i) = [E~ + t, Ef + 1]. □ 

Theorem 7.2. Consider two disjoint balls 1$l(x), B^(y) and £/ie random variables 

:= iBi^r 1 £ V(z;w), TfeH := V(«; w) - z e B L (x) 
zeB L (x) 

(the sample average and fluctuations of the potential in T$l(x)). Let § x be the 
sigma-algebra generated by the random variables {r] y , y £ B^(y); V(z; •), z (jL Bi(x)}. 
Consider the conditional probability distribution function 

^,(*|3 f ») = P{&<«|ffx} 

and zis continuity modulus 

v t„{ s I &c) = SU P ess sup (F ?x (i + s\$ x ) - F£ x (t\!$ x )) ■ 

Suppose that, for some C, C, A, A' , B, B' G (0, +oo) 

Vse[0,l] P{^(s|^x) >CL A s B ) <C'L A 's B '. (7.2) 

P { 3 E e I : xmn{M x (E), M y (E)} > a } < N 2 h(ib) (7.3) 

where 

h{s) :=CL A s B + C'L A 's B ' . 

Proof. Using the decomposition V(z;oj) = £, x {uj) l+ry z (w) in the ball , con- 

sider the respective operator decomposition 

= A x (oj) + £ x (w) 1, A(u) = H + r).(w), 

where, conditional on $ x , the operator A(u>) is non-random. Fix a number b > 
and consider the events Bb{x) (relative to the operator H Bl ^) and, respectively, 
Bb(y); let B = Bb{x) U For any uj $ B, the energies E where M X (E) > a 

are covered by intervals Ji. x = [E~ x , E^ x ], with 5^|Ji |X | < °: and, respectively, 
the energies E where M y (E) > a are covered by intervals J^ y = [E~ , Ef ], also 
obeying J^. |J^ V | < 6. Conditional on $ Xl intervals Ji_ y become non-random, while 
for the intervals Ji, x we can write, by virtue of assertion (B) of Theorem 1 7. 11 

J ilV ( U ) = [l^M E^i^+U")} 

where E^, (u>) are ^-measurable, i.e., non-random under the conditioning by $ x . 
Further, set ti. x = \Ji. x \ and ej. y = \Jj, y \, then 

W ■ Ji,x n Jj tV ± 0} c {\E ( rf\oj) - E^-'°\uj)\ < a, x + e hy ) 

= ||^(w) - < ti,x + ej,j/| 

with -E(w) non-random under the conditioning. Apply the assumption (|7.2[) : 
P{ \£ x (u)) - E(ui)\ < e^ + e^y } <P{ -E{cj)\ <4&} 

< P { v u (4b | &,) > CL A (4b) B } + CL A (Ab) B 
= h(4b). 
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Taking the sum over we obtain the erquired bound: 

P{3E G / : mm{M x (E),M y (E)} > a} < ^P{u; : J h* n J i,i/ ^ i 

< iV 2 /i(4&). □ 

In particular, taking into account Theorem 15.21 we can set, for L = Lk, 

a = a(L k ) = e -T^*) £ *, = L~^\ b = 6(Z fe ) = L^ (1+e)fc . 
These settings give rise to the following corollary of Theorem 17.21 

Theorem 7.3. If there is an integer Lq > 1 and numbers m > 1, a G (1,2) suc/i 
i/iai (|2.7p is fulfilled and 

mm{P , Oo} < CJ 2 Lq K , « > 

i/ien /or some > and a?Z fc > 0, for any interval / Cl with \I\ < 1, 

P{£el: Bi(x) andB L (y) are (£J,m)-5} < 

7.2. On the validity of the assumption (|7.2p . First of all, recall that, by an 
elementary result on Gaussian distributions, if V : Z x 57 —> M is an IID Gaussian 
field, say, with zero mean and unit variance, the sample average £, x of the sample 
{V(z;uj),z G Bl(x)} is independent of the sigma-algebra generated by the "fluc- 
tuations" rj z (uj); moreover, it has Gaussian distribution A/"(0, |Bl(x)|) and admits 
a probability density with Hp^Hoo 

< |B L (x)| 1/2 . In this particular case, Eqn 
(|7.2I) can be replaced by a stronger, deterministic bound: the conditional continu- 
ity modulus (s |3x) is actually independent of the condition and is bounded by 

IbeJU • s. 

Such a situation is rather exceptional, as shows the example of two IID random 
variables Vi(ui), ^(w) with uniform distribution Unif([0,l]). Indeed, in this case 
£ := (Vi+V2)/2, -q — (Vi — V2)/2 and the distribution of £ conditional on rj is uniform 
on the interval of length 0(1— \n\), hence, with constant density 0(|l — 1?7|| ), for 
|?/| < 1; for n — ±1, this distribution is concentrated on a single point. However, this 
example shows also how such a difficulty can be bypassed: excessively "singular" 
conditional distrubutions of the sample mean £ occur only for a set of conditions 
having a small probability. Using this simple idea, Gaume |18j , in the framework of 
his PhD project, established the property (|7.2p for IID random fields with piecewise 
constant marginal probability density. By standard approximation arguments, it can 
be easily extended to piecewise Lipshitz (or Holder) continuous densities, which is 
sufficient for most physically relevant applications. We believe that some variant of 
the property (I7.2[) . perhaps weaker but still sufficient for the purposes of the MSA, 
holds true in a larger class of IID random fields. 

7.3. Exponential spectral localization. The assertion of Theorem 17.31 has the 

same form as in the conventional MSA bound going back to the work by von Drcifus 
and Klein jTJ] (actually, even slightly stronger); therefore, the same argument as 
in [TJ] (having its roots in 17 ) applies and proves that with probability one, all 
polynomially bounded solutions to the equation H(uj)ip = Eip are in fact decaying 
exponentially fast at infinity, thus the operator H{uj) has a.s. pure point spectrum. 
The latter property requires a Shnol-Simon type result on spectrally a.e. polynomial 
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boundedness of generalized eigenf unctions; it will follow by RAGE theorems from 
the dynamical localization proven in Section |SJ 



8. From MSA to strong dynamical localization 

The first rigorous derivations of the dynamical localization from MSA-type prob- 
abilistic bounds on the Green functions have been obtained by Germinet-De Bievre 
[T§] and Damanik-Stollmann 13 . The latter paper had a very eloquent title: " Mul- 
tiscale analysis implies strong dynamical localization" . One of the main ingredient 
of these two works is the analysis of the so-called centers of localization of square- 
summable eigenfunctions; this notion appeared earlier in the work which proved 
instrumental for a number of subsequent researches. Later, Germinet and Klein |20] 
discovered a substantially shorter argument, using more efficiently Hilbert-Schmidt 
norm estimates for spectral projections in a infinitely extended configuration space. 
Formally, [20] considers operators in a Euclidean space R d , but an adaptation to a 
finite-dimensional lattice or, more generally, to a countable graph with polynomially 
growing balls, is quite straightforward. 

In the present paper, working with finite graphs, we bypass the 'hard' analysis 
of spectral projections and replace it by Bessel's inequality. 

The main result of this section can be summarized in the following meta-theorem, 
expressing the 'soft' argument by Germinet-Klein (viz., the finite- volume version 
thereof): " The MSA bounds are essentially equivalent to the strong dynamical lo- 
calization", with a meta-proof: "Apply Bessel's inequality" . 

Owing to the results of Sections EH3 it actually suffices to perform only the 
fixed-energy MSA, even in its simplest form proposed in |26j . 

The extension to an infinite configuration space also admits a short meta-proof, 
going back to the works by Aizenman et al.\ "Apply Fatou lemma" . 

8.1. EF correlators in finite balls. Given an interval J CR, denote by 
the set of all Borel functions <fi : K — > C with supp <fi C I and \\<f>\\oo < E 

Theorem 8.1. Fix an integer L 6 N* and assume that the following bound holds 
for any pair of disjoint balls B^(x), Bl(t/) and some quantity £(£) > 0: 

¥{3E e I : B L (x) andB L (y) are (E,m)-S} < ((L). 

Then for any x,y G Z with d(x,y) > 2L + 1, any finite connected subgraph (of Z) 
Q 3Bj,(i)U Bi(y) and any Borel function <f> € 

E[\(l x \^(Hg(cj))\l y )\] <4e- mL + t(L). (8.1) 

Proof. Fix points x,y £ Z with d(x, y) > 2L + 1 and a graph (/ D Bi(i)UBi()/). 
The operator Hg{uj) has a finite orthonormal eigenbasis {ipi} with respective eigen- 
values {Ai}. Let S = 8Bl{x) U <9Bl(j/) (recall: this is a set of pairs (u,u')); 
note that |S| < 2CjL d , by (|2.1I) . Suppose that for some to, for each i there 
is Zj e {x, y} such that Bi(z;) is (Ai,m)-NS; let {vi} = {x, y} \ {zi}. Denote 
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Hx, y (<l>) = | (la; \^>(Hg(uj))\ ly)\, with fi x ^ y ((j)) < 1. Then by the GRI for the eigen- 
functions, 

fJ>x, v ( < t>) < IMU ^2 \^i{x)i>i{y)\ < ^2 \i>i(Zi)lpi(Vi)\ 

<Y J mvi)\^- mL {ciL d )- 1 

A.ei (u,u')eaB t (ai) 
< e~ mL ]T (C 2 d L d y l I^(«)KI^)I+I^(»)I) 

XiEl (u,u')eS 



using Bessel's inequality and |S| < 2CjL d ) 

< e- mL 2 max (2|| l u || 2 + || l x \\ 2 + \\ l y || 2 ) = ^ mL . 

Denoted = {3E el : B L {x) and B L (y) are (E,m)-S}, withP{S L } < ({L), by 
assumption. Now we conclude: 

E[ Mx , J/ (</))]=E[l SiMx)y (0)]+E[l S c /i:CiJ ,(0)] <C(£)+4e- mi . □ 

It is clear that the exponential term e~ mL can compete with C(L) only in appli- 
cations to the FMM, for the MSA bounds on ((L) are at best sub-exponential in 
L. Otherwise, (,{L) is the dominant term. 

Note also that the decay rate of the bound C,{L) sets natural restrictions on the 
class of the graphs Z, due to the presence of the 'surface' factor |S|, S = S(L). In 
particular, only the FMM-based bounds have the chance to be efficient on trees and 
other graphs with exponentially growing balls. 

8.2. Dynamical localization on the entire graph. Now one can make use of 
a simple argument employed earlier by Aizenman et al. [HO], in the framework of 

( H\ 

the FMM which always starts as a fixed-energy analysis. The quantities /ix,y (4>) = 
(l x | 4>(H) | l y ) defined, for example, for bounded continuous or Borel functions cj), 
generate signed (i.e., not necessarily positive) spectral measures associated with a 
self-adjoint operator H : 



J dn x H y \E)${E) := (1,10(1011, 



In particular, we can consider, with x,y,u G Z fixed, measures fi x related to 
operators H^ L (u), for all k > 0, as well as their counterparts fj, XtV for the oper- 
ator H on the entire graph Z. A sufficient condition for the vague convergence 
fi x y — > fj, x>y as k — > oo is the strong resolvent convergence H Bl ^ — > H. Such 
convergence is well-known to occur for a very large class of operators, including 
(unbounded) Schrodinger operators in Euclidean spaces and their analogs on the 
so-called quantum graphs. Indeed, for (not necessarily bounded) operators H n with 
a common core V to converge to an operator H with the same core, it suffices that 
Hk^P Hip strongly for any element ijj £ V (cf. [35]). For finite- volume operators, 
one can usually find an appropriate core T> formed by compactly supported func- 
tions ijj; for finite-difference Hamiltonians on graphs (even unbounded, e.g., for DSO 
with unbounded potentials) one can choose as T> the subset of all functions with 
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finite supports. On such functions, H& L („)^/> — > Hip as k — > oo (by stabilization), 

therefore, the spectral measures converge vaguely: /i* — > /i x ,j/- By Fatou lemma, 
for any bounded Borel set icl, one has 

\t*zM\ <liminf \t£, y (A)\ 

(here := sup{/x(0), |j0|| < 1, supp (f> C A}). Taking the expectation and 

using the uniform upper bounds on EF correlators in finite balls, we conclude that 



E 



sup \{l x \((>(H(u))\l t 



(using the inequality £~ K *- 1+6 '' > < Ce~ alnl+C Lk , C,a,c > 0). In particular, with 
functions <f>t '■ A i-» e~ ltx , we obtain the strong dynamical localization property for 
the ensemble of random Hamiltonians H(lu). 

9. Sub-exponential bounds on EF correlators without bootstrap 

Now we will show how the polynomial (or slightly stronger than polynomial) 
decay bounds from Section [5] can be substantially improved and replaced by sub- 
exponential ones. Germinet and Klein proved in Q0' a highly optimized and very 
general sub-exponential decay bound for a large class of random differential oper- 
ators (an adaptation to lattices and graphs is straightforward). Unlike [3D], our 
aim here is to obtain an elementary proof in the simplest situation, without a more 
involved bootstrap procedure. In the light of Sections [6HH1 it suffices to work with 
the resolvents at a fixed energy E E I C M.. 

The advantage of the method presented below is that it gives rise to exponen- 
tial bounds on the decay of eigenf unctions, while using an exponential sequence of 
scales, ife+i = YLk, as in j2S] and in [20] . gives directly only a sub-exponential 
bound. (Recall that [20] uses several multi-scale analyses to obtain final results, in- 
cluding exponential spectral localization. In |26j . it was indicated that exponential 
localization requires scales Lk ~ Lq .) 

The main idea of the method described below is quite natural. The MSA in- 
duction shows clearly that the exponent k > in the power-law bound of the form 
P { B^ k (x) is (E, m)-S } < L~^ K grows with the number Kk of allowed singular Lfc-i- 
balls inside B^ k (x). We allow for a growing number Kk ~ L c k , c € (0, 1), and use an 
elementary probabilistic bound on such an event, close in spirit to the Poisson limit 
theorem. In |20j . a similar effect is achieved by a refinement of an idea from |26j : 
replacing the sequence of scales Lk ~ Lq , a > 1, by a slower growing sequence 
Lfc ~ Y k Lo, while keeping uniformly bounded the maximal number Kk of allowed 
singular cubes. 

9.1. Multiple singular balls: a probabilistic estimate. 
Lemma 9.1. Suppose that for any ball B^ j (a;) C B^ j+1 (u) one has 

P { B Lj {x) is (E, m)-S } < e- L i , 5 > 0. 

LetM(uj) be the maximal cardinality of collections of pairwise disjoint (E,m)-S balls 
C = {B^^Xi), i = 1, . . . , Af}. Then, for a > 6 and Lq is large enough, 
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Proof. Fix a possible (unordered) configuration of centers Xi of disjoint (E, m)- 
S balls, i — l,...,k. Let N — |Bl. +1 (u)|. The number of such configurations 
for a fixed k is bounded by N(N - 1) • • • (N - k + l)/k\ < N k /k\, since choos- 
ing every center in the sequence x\,X2, ■ ■ ., excludes at least one possible position 
for the next center (indeed, many more). For a given configuration, the events 

{Bl (xi) is (E,m)-S} are independent, with probabilities < p :— e Lj , so 

k=n k=0 v ' 

for p < N- 1 and n > 3. With JV < CLf 1 and n := [Lf a ~ l \ one has TVp < 

P{ > n } < (Np) n < exp {- (i| - ClnL^) (if"" 1 ' - l) } 
<exp|--i i+1 «" 

The condition s+a ( a ~ 1 ) > <5 is equivalent to the assumed inequality a > 5. There- 
fore, for some 5' > 5 and Lq large enough 

P{ N > Lf a - 1] } < e-i £ '+i < ^e" L i +1 . □ 

9.2. Decay of (£, q)-subharmonic functions with "singular" points. The ra- 
dial descent bound given by Lemma 14.31 will require an adaptation. 



Definition 9.1. Let Q be a finite connected graph and L > £> two integers and 
q £ (0, 1). Consider a ball Bl{u) C Q and a function / : Q — s- M + . 

(1) We say that a point i £ Bi(«) is (£, q)-regular for the function / iff 

f(x) <q max f(y) 

and denote by 1Z(f)(c Bl(u)) the set of all regular points of function /. 

(2) Given a point x e Bj,(u), let R(x) be the smallest integer such that S_r( x ) C 
7£(/); if no such integer exists, we set formally R(x) = +oo. 

(3) We say that / is (£, q, 7\L)-subharmonic in Bl(u), with 1Z — TZ(f), if for any 
point x with R(x) < oo and for all r > such that 

S r+e+ i(u) C B L+ i(u), S r (u)c7£, 

one has 

/(a:)<g max /(„). (9.1) 

yeB r+e+1 (x) 

Lemma 9.2. Let a function f : Q — > R + , defined in a finite connected graph Q , 
be (£,q,1Z)-subharmonic in Bi(ti) C Q. Suppose that the set 1Z C = Bl(u) \ 1Z is 
covered by a family of annuli 

A = {Ai,l <i<n}, A, = B 64 (u)\B 0i (u), h - Oj < Ci£,a € N*, 

of total width 

w(A) = ^(6, - o<) < c ^ = C ^ 



20 



V. CHULAEVSKY 



ci = 2 

Ji = [ax,h] h(r)<qh(r>) 



h — i — ^ ^ — i — i — i — • — ♦ — i — h 



(2 + l)(£ + l) 

<2. 2 (£ + i) 




Figure 2 . Example for the proof of Lemma 19.21 Recursion for 
r € 7Z (step of length £ + 1) and for s G 7?. c (step of length < | J2I). 



wrf/i CeN and 2C(^ + 1) < L.Then 

f(u)<qim\~ 2C M(f,G). (9.2) 



L + l 



intervals of the 



Proof. Divide the integer interval [0, L + 1] n N into N := 
form Ij = [j(£ + l), (j + 1)(£+1) — = 0, 1, n— 1, and, eventually, a remainder 
which will be unused in the argument. Call an interval Ij good if S^ + i) C 1Z. 
Since B^(m) C £7, the sphere Sl+i(u) is non-empty. The radial projection 

Bl+i(m) 3i4 d(u, x) € [0, L + 1] 

maps an annulus of width Ci£ onto an interval of length Ci£, covered by at most 
+ 1 < 2ci adjacent intervals of the form Ij. Therefore, the entire set TZ C is 
radially mapped onto a subset of [0, L] covered by a family of at most 2C intervals 
Ij . Respectively, at least K = N — 2C intervals Ij 1 ,Ij 2 ,..., Ij K must be good, and it 
follows from the hypotheses that K > 1. Further, let h(r) := max ie B r (u) f(%) > 0, 
r G [0,L + 1]. Note that if x £ Bj(^ +1 )(u) and z € B^ + i(x), then z £ B( J+1 )^ +1 (u). 
For x € with d(w, x) =: r € one can apply (|4.1D . so for any 1 < i < K, one 

has 

h(ji(£ + 1)) = max f(x)<q max /(z) 
a:eB Ji(f+i)( u ) zeB («+i)(«+i)(™) 

= ^((j 4 + l)(£ + l)). 
The value i = K, i.e., the index jx, is admissible here, since 



(9.3) 



(j K + l){£ + 1) < if (£ + 1) + 2C(^ + 1) < X + 1, 

and the set Sl+i(u) is nonempty, by assumption. For i — K, we have, by mono- 
tonicity of the function h, the following upper bound on the RHS of ([9. 3ft : 



qh((ji + l){£ + 1)) < qh(L + 1) < gH/]^, 
so by backward recursion i = K — 1, 1, one obtains in if — 1 steps: 

/(«) = Mo) < h(h(£ + i)) < q K ll/iu < L^J- 2C ii/Hoo. □ 

Definition 19.11 is tailored so as to suit the following analog of Lemma 14.21 
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Lemma 9.3. Consider a finite graph Q, a DSO Hg and a ball B^ k +1 (u) C Q. 
Fix some E € K, assume that Bi fc+1 (u) is E-CNR, and let 1Z = 1Z(G,E) be the 
set of points x £ Q such that the ball Bi, k (x) C Q is (E,m)-NS. Then for any 
V £ G \Bi fc+1 (u) 7 the function x i-> \G-B Lk+i (u)( x ,y, E)\ is (Lk,q,7Z)-subharmonic 

in Bi, k+X (u), iirat/i q = e m< - 1+ ^ L fe ) Lfc 7 */£o *s ^r9 e enough. 

Proof, li x £ 7Z, then, by the m)-NS property of ~BL k (x), 

f(x)<e-^ m < L ^ L »M(f,B Lk+1 (x)). 

For a; £ 1Z C with < oo (the radius R(x) is defined as in Definition 19. ip . we 

have, by E-NK property of the ball Rr( x )-i{u), stemming from the assumed E'-CNR 
property of B Lk+1 (u), 

f(x)< \8B Lk (x)\ WG* m(x,v;E)\\ „ max,, JG Bi _ M (v',y;E)\ 



d(u,v')=R(x) ' TBi fc + i (u)l 



\G 



BL k + 1 (u) 



(v',y;E)\ 



< e k max 

d(u,v')=R(x) 

and since SW,,.) (w) C TZ, 

<e L " max e^ 7(m ' ifc)Lfc X(/,B 1+Lfc (w')) 

< e-^ m ' L ^ +L "M(f, B R(x)+Lk+1 (u)). 

For L Q large enough and r < 1 — /3, one has 7(m, Lk)Lk — EP k > m(l + jL' k T )Lk. □ 

9.3. Scaling with sub-exponential bounds. From this point on, we fix the key 

parameters used in the scale induction as follows: 





= * 


5 = i < /3 


cr = | > <5 







Definition 9.2. A ball Bl h1 (i) is called -©-tunneling (E-T) if, for some E £ I, it 

pairwise disjoint (-E 1 , m)-S balls of radius 



contains a collection of more than E 



r(a-l) 



ifc, and iS-non-tunneling (iS-NT), otherwise. 

Lemma 9.4. // a 6aZZ B L(c+1 (w) is E-NT and E-CNR, then it is (E, m)-NS. 



Proof. Fix two points x,y £ B^ k+1 (u) with d(x, y) > L^ +g and let R = d(x, y) — L 



By Lemma 19.31 the function / : z >->• G-Q Lk+1 {x, z) is (Lfe, q, 5)-subharmonic in 
B = Br(x), with q = e~ mi - 1+ i L k ) ifc , and the E-NT assumption guarantees that S 
can be covered by at most L a ^ a ^ balls of radius 2Lk, hence by a collection A of 
annuli centered at u of total width w(A) < L^ a l ' ■ 4L&. Therefore, 



R-u,-(A)-2L k 

f(x) < q W H/ll 
Recall that g — a{a — 1) = 2r, so 



R-5L 



1+ct(qi-1) 



>Rll-5L 



< 



~ Il/I 

= R(l - bL 



(9.4) 
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Put this lower bound into (19.41) and write, with the convention — InO = +00: 



(with m > 1, i? = d(x, y) — Lk > L^ e — Lk > \L X k +B and Lq large enough) 

> mR ((1 + (1 - 6^ 2T ) - 2L- k 1 -^) 
(use /? = \, a/3 = | => 1 + p - a/3 = § > 2r) 

> mR (1 + §LjT ) (1 - 7L^) > m (l + ±^7 r ) d(z, t/) 
>7(m,ife+i)d(a;,y). □ 

Consider the following property which we will prove by induction for all fceN: 
S(k,E): For any ball Bi fc (x) C Z, one has 

P { B Lk (x) is (£, m)-S } < e _i * . (9.5) 
Theorem 9.5. S(k,E) implies S(k + 1, E) . 

Proof. Denote by N(Br lk+1 (x)) the maximal cardinality of collections of pairwise 
disjoint (E,m)-S balls of radius Lk inside Br, k+1 (x). Introduce the events 

B k+1 ^{N(B Lk+1 (x))>Ll {a - 1) ], 
£ k+1 = { B Lk+1 {x) is E-PR}. 

By Lemma [Hill {w : Bi fc+1 (x) is (E 7 m)-S} C Sk+i U Bk+i- By Wegner estimate 
(|2.8p . using (i > 5, we have 

P{^+i}<|e-^<ie-^ 

so the claim follows from Lemma saving that P{6 fe+ i} < ie- L ^. □ 

For the initial scale bound S(0) (indeed, any desired probabilistic bound, cf. jl2j 
to hold true, it suffices to pick \g\ large enough. Therefore, we come by induction 
to the following 

Theorem 9.6. Assume that the random potential fulfills the regularity conditio-^ 
(|2.7[) . Then for all \g\ large enough, the property S(k) holds true for all k > 0. 

Remark 9.1. A more tedious (but elementary) parametric analysis shows that one 
can actually choose the exponent S £ (0, 1) arbitrarily close to 1, thus getting a 



2 As was said, the Lipshitz continuity condition (12.71) can be relaxed to the Holder continuity. 
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sub-exponential decay very close to the exponential one. Indeed, the complete set 
of requirements for the scaling parameters is as follows: 

< p < a - 1 

< a < — - — , < 6 < min{0, <t} 
a — 1 

< 2r = min{p — a(a — 1), 1 + p — a/3} 
A direct inspection shows that for any e € (0, j), one can set, e.g., 



a = 1 + 4e 


p = l-e 


6 = 1 - 2e < (3 


a = 1 - e > 8 


p = 4e- 2e 2 


T = e 2 = \{p-a{a-l)) 



10. ON A DIRECT APPROACH TO THE DECAY OF EIGENFUNCTIONS 

We have seen that the fixed-energy analysis of resolvents implies the spectral and 
dynamical localization. On the other hand, it is also possible to adapt the approach 
from the Spencer's work |26j to a direct analysis of eigenfunctions in arbitrarily 
large finite balls. For single-particle models, such an adaptation has been proposed 
in our recent paper [8]. The new scaling scheme can be summarized as follows: 

• The key notion becomes that of an m-localized ball. We say that a ball 
Bl(u) is m-localized if the eigenbasis {ipi} of the operator H Bl ^ fulfills 
the following condition: for all points x : y e Bi(«) with d(x,y) > L 7 / 8 , 

£ Mx)My)\<^ (m ' LWx ' v} - 

A*GE(J? BlW ) 

• A ball B^ fc+1 (u) is called E'-completely non-resonant (_E-CNR) if it is E'-NR 
and contains no E-K ball of radius > L\. . 

• It is readily seen from the eigenfunction expansion of the resolvent that an 
m-localized ball which is E'-NR must be (E,m)-NS. 

• A direct analog of Lemma 15.11 is still valid. 

• Consider the bounds (which we will denote by LOC(fc), k > 1) of the form 

P{B Lk (u) is m-localized } > 1 - L^ K(1+fl) *. (10.1) 

The initial scale bound LOC(l), for any m > 1, Lq > 1 and \g\ large 
enough, is easily obtained by elementary perturbation theory for self-adjoint 
operators (which are finite-dimensional here) with simple spectrum. The 
induction step (LOC(fc) =>• LOC(fc + 1)) is performed as follows. 

— Assume LOC(fc) and consider a ball B£ fc+1 (it). If it is not m-localized, 
then by an analog of Lemma 15.11 for some E € M it must contain two 
disjoint (E,m)-S balls Bi, k (x), Bi fc (y). One can easily infer from the 

_r/3/2 

Wegner estimate that with probability > 1 — e k + 1 , for any E e K 
either B Lk (x) or B Lk (y) is E-CNR. 

— If BL k (x) is (E,m)-S and i?-CNR, it must contain a pair of disjoint 
m-non-localized balls of radius Lk—i- By virtue of the inductive as- 
sumption LOC(fc), the probability of the latter event is bounded by 

n r 2da 2 -2n(l+e) k - 1 
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— For a e (1, V2), k > I" a and an appropriately chosen 6 > 0, the 
above bounds imply LOC(fc + 1). 
The property LOC(/c), proven for all k > 0, is already a form of localization 
of eigenf unctions. In addition, it implies the usual variable-energy MSA bounds, 
hence, the strong dynamical localization and an exponential decay of eigenfunctions. 

It is to be emphasized that the main analytic tool of the simplified MSA remains 
the elementary Lemma 14.11 (and Lemma 14.31 easily stemming from it). As was 
mentioned earlier, the idea of the "two-sided" estimates of the Green functions 
(and, similarly, eigenfunctions) goes back to Spencer's work |26) . 

Conclusion 

We have shown that the fixed-energy probabilistic analysis of random Anderson- 
type Hamiltonians implies, in a fairly general and elementary way, stronger manifes- 
tations of the Anderson localization phenomenon, viz.: spectral and strong dynam- 
ical localization. The new method, going in the same direction as the well-known 
Simon- Wolf criterion of localization, applies also to multi-particle systems (as will 
be shown in our forthcoming work) for which no analog of the Simon- Wolf approach 
has been developed so far. Moreover, combined with the simplified Germinet-Klein 
argument, it gives rise to the strong dynamical localization, not only spectral local- 
ization. Therefore, a very simple scaling procedure going back to Spencer's work 
[2"5] . as well as its counterpart for interacting mult i- particle systems, results in a 
simple proof of strong dynamical localization for a large class of models. 
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